Dispersive contributions to e + p/e p cross section ratio in forward regime 

Mikhail Gorchtein 
California Institute of Technology, Pasadena, CA 91125, USA 
(Dated: February 2, 2008) 

Two-photon exchange (TPE) contributions to elastic electron-proton scattering in the forward 
regime are considered. The imaginary part of TPE amplitude in these kinematics is related to the 
DIS nucleon structure functions. The real part of the TPE amplitude is obtained from the imaginary 
part by means of dispersion relations. We demonstrate that the dispersion integrals for the relevant 
elastic ep-scattering amplitude converge and do not need subtraction. This allows us to make clean 
prediction for the real part of the TPE amplitude at forward angles. We furthermore compare e + p 
and e~p cross sections which depends on the real part of TPE amplitude, and predict the positron 
cross section to exceed the electron one by a few per cent, with the difference ranging from 1.4% 
to 2.8% for electron lab energies in the range from 3 to 45 GeV. We furthermore predict that the 
absolute value of this asymmetry grows with energy, which makes it promising for experimental 
tests. 
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I. INTRODUCTION 

Recently, much attention has been attracted by the 
two-photon exchange (TPE) contribution to the elastic 
electron-proton scattering. On the one hand, the recent 
experimental data on Ge to Gm ratio at higher momen- 
tum transfers show significant discrepancy between the 
results obtained with the new polarization transfer tech- 
nique 0] and those obtained using Rosenbluth separation 
0- Since the TPE contribution is the only largely un- 
known order a em correction to the elastic ep-scattering, 
it was argued that a proper inclusion of these effects into 
the analysis of the ep-scattering observables would rec- 
oncile the two measurements |3|. Recent theoretical cal- 
culations of the real part of the TPE amplitude in dif- 
ferent models seem to support this idea, though at the 
qualitative level Q|, p|, |g- At the moment, it is only the 
clastic (nucleon) intermediate state contribution together 
with the associated with it IR divergent part, needed 
for an analysis of the experimental data, that one can 
be confident about PI. The partonic model "handbag" 
calculation of Ref. though definitely represents im- 
portant physics contribution at high energy and momen- 
tum transfer, suffers of unphysical IR divergencies which 
can be presumably cancelled by including the "cat ears" 
mechanism. As for the further inelastic contributions, 
only estimates with the A(1232) for the real part of the 
TPE amplitude exist in the literature 

Another experimental test of TPE effects is proposed 
at JLab and at VEPP-3 storage ring @, where the 
e + p/e _ p cross section ratio will be studied in a wide 
angular range. A deviation from 1 originates to lead- 
ing order in a em from the interference between the one 
photon exchange contribution which is linear in lepton 
charge, and the real part of the TPE contribution, which 
is quadratic in lepton charge. 

On the other hand, over past 5 years new observables 
become accessible experimentally, the single spin asym- 
metries with beam electron or target (recoil) proton po- 



larized normally to the reaction plane [9( . These observ- 
ables have been first studied theoretically back in 1970's 
[Tof and were shown to be directly sensitive to the imag- 
inary part of the TPE amplitude. Furthermore, all IR 
divergent terms cancel in these asymmetries (unlike the 
cross section), which makes them especially attractive 
since they offer the most clean experimental tests of the 
TPE effects. There is a considerable interest to this class 
of observables from the theory point of view, as well [Tl| . 

The TPE amplitude is related to Compton scattering 
with two space-like photons (virtual-to-virtual Compton 
scattering - VVCS) which is largely unknown. For the 
inelastic intermediate states contributions, it is only the 
imaginary part of the VVCS amplitude in the forward 
direction that is well studied, since it can be expressed in 
terms of the DIS structure functions. So, the only kine- 
matical point where a firm prediction for the TPE am- 
plitude can be made is the exact forward limit. Unfortu- 
nately, all the observables which measure the TPE effects 
necessarily vanish at that point. It was proposed 
to combine the DIS input in forward direction with the 
phenomenological t-dependence taken from that of the 
Compton scattering differential cross section which 
is measured at high energies and low values of t. Mod- 
elling in this way the imaginary parts of the electron 
helicity-flipping ep-scattering amplitudes, it was possible 
to successfully describe the data on beam normal spin 
asymmetry at forward angles [l3T |. The motivation of 
this work is to extend the phenomenological approach 
of the previous paper and compute the imaginary parts 
of the electron hclicity-conserving amplitudes. The real 
parts which enter the expression for the elastic cross sec- 
tion and the charge asymmetry are obtained through dis- 
persion relations. Combined with the phenomenological 
^-dependence which proved itself plausible for the imag- 
inary part, it will allow us to predict the high energy 
dispersive contributions to the e + p/e~p cross section ra- 
tio and the elactic cross section. 

For the case of low electron energies, a dispersive ap- 
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proach was adopted in The authors of that work 

used once-subtracted dispersion relations in the annihi- 
lation channel t at fixed energy « 0. At this energy, the 
hadronic part of the TPE graph may be parametrized 
through the nucleon polarizabilities. To calculate the 
subtraction constant, the authors notice that the imagi- 
nary part of the hadronic amplitude at t — is given in 
terms of DIS structure functions. They furthermore use 
dispersion relations for these structure functions, and by 
doing the low energy expansion express the correspond- 
ing real parts as moments of the structure functions. Af- 
ter that, these real parts are embedded into the TPE 
amplitude. 



Instead, our approach is designed for high electron en- 
ergies, therefore no connection with the polarizabilities 
can be made, neither the low energy expansion can be 
performed to express the final result as an integral over 
the DIS structure functions' momenta. First, we calcu- 
late the imaginary part of the TPE amplitude in terms 
of the electron phase space integral over DIS structure 
functions and obtain the imaginary parts of the invari- 
ant ep-scattering amplitudes. Instead of using dispersion 
relations for the hadronic amplitude separately, we pro- 
ceed with the dispersion relations for the invariant am- 
plitudes for elastic ep-scattering and demonstrate that 
they converge and no subtraction is needed. In this way, 
we obtain a clean prediction for the real parts of these 
amplitudes and the observables which they enter. 



The imaginary parts of helicity-flip amplitudes were 
shown to be dominated by the region of low virtualities 
of the exchanged photons, and it prooved adequate to 
approximate the DIS structure function W\ by its value 
at the real photon point, where it is related to the total 
photoabsorption cross section. However, this approach is 
not expected to work in the case of the helicity-conserving 
amplitudes, especially for their real parts where all the 
values of these virtualities are allowed. Therefore, we 
include in our calculation the phenomenological form of 
the Q 2 dependence of the virtual photon cross section ob- 
tained from low- a; data. For the moment, we will concen- 
trate on the observables averaged over the nucleon spin, 
i.e. the cross section and the beam normal spin asymme- 
try. We leave the consideration of the target normal spin 
asymmetry to an upcoming work. 



The paper is organized as follows. In Section [HI we 
overview the elastic ep-scattering amplitude and the kine- 
matics of the reaction. The calculation of the imaginary 
parts of the invariant ep-scattering amplitudes is per- 
formed in Section IIIII These imaginary parts are then 
used as input for the dispersion relations to determine 
their real parts in Section IIVI and the convergence of 
these dispersion relations is discussed. In Section we 
present the results of the calculation. 



II. ELASTIC ep-SCATTERING AMPLITUDE 

In this work, we consider elastic electron-proton scat- 
tering process e(k) + p(p) — > e(fc') + p(p') for which we 
define: 



P = 



K = 



P + P 
2 

k + k' 



q = k — k! — p — p, 



(1) 



and choose the invariants t = q 2 < 1 and v = (P-K)/M 
as the independent variables. M denotes the nucleon 
mass. They are related to the Mandelstam variables s — 
(p + k) 2 and u = (p — k') 2 through s — u — AMv and 
s + u + t = 2M 2 . For convenience, we also introduce 
the usual polarization parameter e of the virtual photon, 
which can be related to the invariants v and t (neglecting 
the electron mass m): 



e = 



u 2 ~M 2 t{1+t) 
v 2 + M 2 t{1 + tY 



(2) 



with r = — t/(4M 2 ). Elastic scattering of two spin 
1/2 particles is described by six independent amplitudes. 
Three of them do not flip the electron helicity |3|], 



T, 



no flip 



— u(fc')7 M u(£;) 



(3) 



F 2 — 

M 



while the other three are electron helicity flipping and 
thus have in general the order of the electron mass m e 
El: 



T 



flip 



m e e 



F G u(k')j 5 u(k) ■ u(p')j 5 u(p) 



(4) 



In the one-photon exchange (OPE) approximation, two 
of the six amplitudes match with the electromagnetic 
form factors, 



F. 



3,4,5,6 



(M) = 



mi 





(•5) 



where Gm (t) and F 2 (t) are the magnetic and Pauli form 
factors, respectively. For further convenience we define 



1 In elastic ep-scattering, the usual notation for the momentum 
transfer is Q 2 = — q 2 but we prefer the more general notation t 
to avoid confusion with the incoming and outgoing photon vir- 
tualities. 
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also Ge — Gm — (1 + t)F 2 and F\ = Gm — F 2 which 
in the Born approximation reduce to Sachs electric form 
factor Ge and Dirac form factor F\, respectively. It is 
useful to separate one- and two-photon exchange effects 
explicitly, 



Gm 
G e 
F 2 



Gm + 5Gm, 
Ge + SGe, 
F 2 +SF 2 , 



(6) 



where Gm, Ge, and F 2 are the usual form factors, while 
the two photon effects are contained in the quantities 
SGm, 5Ge, SF 2 . 

The unpolarized cross section is 



da 



TCTR 



dan 



dflLab £(1 + T) dflLab 

with the usual Rutherford cross section 
da 



4n -C0S 2 f E' 3 

~~E' 



dil 



Lab 



(7) 



(8) 



O the electron Lab scattering angle and E(E') the incom- 
ing (outgoing) electron Lab energy. The reduced cross 
section an is given by 



a R = G 



£ ^2 
T 



2-G E Re [5Ge 



2G M Re (SG M 

0(e 4 



e — F3 
M A 



-F 3 

M 



(9) 



The OPE contributions to the form factors depend lin- 
early on the lepton's charge, while the TPE contributions 
are quadratic in the lepton's charge. Therefore, the in- 
terference terms between the OPE and TPE will enter 
e~p and the e + p cross sections with opposite sign. The 
e + p/e~p cross section ratio is thus given by 



R e+/e- 



a e + 



e~<-p 



a&— p 
1 - 



(10) 



Gm Re (SG M + e^F 3 



-G^Re (sGe 



-h 

M 3 



For a beam polarized normal to the scattering plane, 
one can define a single spin asymmetry, 



the early 1970's ^(j that such asymmetries are directly 
related to the imaginary part of the T-matrix. Since 
the electromagnetic form factors and the one-photon ex- 
change amplitude are purely real, B n obtains its finite 
contribution to leading order in the electromagnetic con- 
stant a em from an interference between the Born ampli- 
tude and the imaginary part of the two-photon exchange 
amplitude. In terms of the amplitudes of Eas. (|3l4|) . the 
beam normal spin asymmetry is given by: 



-^ v / 2^r^)VTT7(rGl f +eG 



tG m Im-Fa + G E lmF 4 



Fi^lmF 5 ] (13) 



III. IMAGINARY PART OF THE FORWARD 
ELASTIC ep-SCATTERING AMPLITUDE 




FIG. 1: Imaginary part of the 27-exchange diagram 

The imaginary part of the diagram in Fig. \T]is given 
by the integral 



ImT; 



2- 



^riW'i^. (14) 



(2^)32^ q\q\ 
where the leptonic tensor is given by 



(15) 



The imaginary part of the spin-averaged part of the 
hadronic tensor is expressed in terms of the DIS structure 
functions W\ and W 2 . We will make use of Callan-Gross 
relation between them, thus we obtain 



cr T + a l 



(11) 



where a^ (ai) denotes the cross sesction for an unpolar- 
ized target and for an electron beam spin parallel (an- 
tiparallel) to the normal polarization vector defined as 



5 n 



0, 



[k x k'] 
\kx k'\ 



(12) 



normalized to (S ■ S) — —1. Similarly, one defines the 
target normal spin asymmetry T n . It has been shown in 



fXV 



2-kWx \ -9' 

(qi ■ 92) 



P»q\ + P v q^ 



PK 



{P-Kf 



(16) 



with q^ = k — k\ the incoming and qt^ = k' — k\ the 
outgoing photon momenta and their average K — 2i±2a 

We next contract the leptonic and hadronic tensors and 
keep the lepton mass non-zero in order to cross-check the 
calculation with the previous work for the beam normal 
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spin asymmetry |l7| . The result reads 



Lu ■ ImlVf"' = 2ttW 1 i m e uu 



2Pfci 



PK 



_, 2K 2 Pki I Pk x PK PK 

+u ^[{pky [JpW + pk 

Ql(2PK + PK_ 2 
PK \ PK PK 



(17) 



In the above equation, we used the notation Q 2 = 
—q\. Throughout the calculation, we neglect 



terms 



tt^ which are small in forward kinc- 



(PK) 2 

matics. 

The above formula represents the TPE amplitude in 
the nucleon helicity-conserving channel, averaged over 
the spin projections. The dependence on the electron 
spins is retained. Our goal is to obtain the expressions 
for the invariant ep-scattering amplitudes defined in Eqs. 
00}. For this, we will need to restore the dependence 
on the nucleon spin. It can be done as follows. There 
exist four independent scalars which can be formed from 
the initial and final nucleon spinors with the 7-matrices 
and the four- vectors that fix the external kinematics. 
They are N'N, N'fSST, N'j 5 f<N and N'j 5 N. 2 However, 
the fourth structure vanishes in the exact forward limit 
and in the helicity conserving channel. Furthermore, the 
third structure drops after averaging over spins (that is, 
it is non-zero in the GDH sum rule-like cross section dif- 
ference r 1 /2 i A';i/2,A-7 1 _ 1 / 2i A';-i/2,A, but not in the sum). 
Finally, we are left with only two structures N'N, N'fiN. 
By a direct calculation, it can be shown that in the for- 
ward limit, 



1 



1 = - 



A I 



2^1 s-M 2 

spins 



N'N - 



s + M 2 

(s - M 2 ) 



; N'fN } (18) 



Therefore, we now modify the hadronic tensor by sub- 
stituting 



M 



M 2 



N'N - 



M 2 



{s - M 2 ) 



: N'KN} (19) 



This procedure allows one to restore the nucleon spin 
dependence of the hadronic amplitude from its spin- 
averaged part in a correct way, as long as only cross 
section channel is used as the phenomenological input. 



2 The vector N'y^N, the axial vector N'^^^N terms, as well 
as the tensorial terms N'a^N and N'a^'yf.N in the adopted 
formalism have to be contracted with the four independent 
four- vectors which form a basis, that is K^,P tlL = P M — 
-^if^q^iV = e^^P^Kpqj. Working out the Dirac al- 
gebra, it can be shown that the only independent structures are 
N'N, V'/iV. N'^fN and ^'75^ 



We can now identify the invariant amplitudes Fi, i 



,6 comparing Eqs. (|3I4I17I19[) : 

ImG M = 



(20) 



ImF 2 = 

ImF 3 = 

ImF 4 = 

ImF 5 = 
lmF 6 



-2ne 2 t 



M 2 



M 2 



s-M 2 



M 2 



ImF, 



2ne 2 t 



M 



s - M 2 
M 2 



M 2 



lmF 4 



= 0, 



where we introduced shorthands 
{W 1 (w 2 ,Q 2 ,Q 2 )} = 
2K 2 Pk 1 ( PhPK 



1 



(PK) 2 \ (PK) 2 

Q 2 ( 2PK PK 
PK \ PK PK 

W 1 (w 2 ,Q 2 ,Q 2 ) = 



(2n) 3 2E l Q\Q\ 
PK 



- 2 



W^^lQj) 

d 3 A?i 1 
(2nf2E 1 Q\Q 2 



2Pfc i W 1 2 



i) 



(21) 



Note that the fact that ImGM vanishes implies that 
ImFi = — Imi^. In the above formula, the first argument 
of W\ is the invariant mass squared of the 7*p system, 
w 2 = (P + K) 2 . In the integral J we have to keep the 
term ~ K 2 = — | till the end, since it is needed to cancel 
the 1/t behaviour of the integral I\ introduced below. 

The next step is to perform the phase space integrals 
in Eq. 121|l . First, we express the nucleon structure 
function W\ through the virtual photon cross section 
a 1 * p (w 2 ,Q 2 ): 



Wi 



w 2 - M 2 
2ire 2 



>7*p 



(w 2 ,< 



(22) 



We will need an explicit form of energy and Q 2 de- 
pendence of this cross section. We will use the phe- 
nomenological form proposed in Ref.JllJ which amounts 
in factorization of the virtual photon cross section into 
purely u> 2 -dependent part, and the low- a; scaling variable 
^-dependent part: 



a Tp (w 2 ,Q 2 ) 



Regge 
IV 



(w 2 ) 



I[ V (A 2 (w 2 )),r ]0 ] 



(23) 



where a^" e (w 2 ) is the total photoabsorption cross sec- 
tion in Regge regime and is conveniently parametrized in 
terms of two Regge trajectories, 



Regge 
IV 



(w 2 ) = A p {w 2 ) a >- 1 + A P (w 2 Y 



(24) 
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with A p = (145±2)/i6, A P = (63.5±0.9)/A a p = 0.5 and 
ap = 1.097± 0.002 the parameters of the p and pomeron 
trajectories, respectively. In the above formula, w should 
be taken in GeV. The scaling variable 77 is defined as 



V 



Q 2 



A 2 (w 2 ) ' 



(25) 



Fur- 



with its minimal value 770 = Tj{Q = 0) = ■ 

thermore, the function of energy A(u> 2 ) is determined by 
a fit to DIS data which are well reproduced by 

k\w 2 ) = atf+wlf; 



(26) 



with Ci = 0.34 ±0.05, C 2 = 0.27 ±0.01, w 2 = 882 ±246 
GeV 2 , and to 2 , = 0.16 ± 0.01 GeV 2 . The general form of 
the function 1(77,770) can be found in |16|. Here, we will 
use its asymptotic form for low values of 77, 



= In ( - 

,11 



(27) 



Furthermore, we divide the virtual photon cross section 
in two parts, the Q 2 independent and the rest, 



Regge 
IP 



(w 2 )) 



IP 



Regge 
IP 



lnfl 



Q 2 /m 2 ) 



In 



(K) 

\ m oJ 



(28) 



A. Electron helicity-flip amplitudes 



We start with the integral J B . Introducing the di- 
mensionsless variable z = ^ and consistently neglecting 
terms ~ t/s, we can write the integral J g as 



(s - M 2 )E 2 
167r 4 e 2 



(29) 



\ / dzz 2 {\- z)a 7P (M 2 + {s- M 2 )(l- z)) 



J QiQl 



i ln(l + Q 2 /m 2 ) 

bx(K) 



The integral corresponding to the first term in the 
square brakets was performed before (see H3)E3)- ^ n or " 
der to obtain the correct i-dependence, we keep the pho- 
ton indices in the denominator: then, the first term in the 
square braket leads to the leading 1/t behaviour. For the 
second term, the denominator can be written as l/(<5 2 ), 
as this corresponds to neglecting the order t 2 corrections 
to the leading term. The second term can be related to 



the table integral / 



111(1+3;) 



dx 



lnx 



l+x 



ln(l 



The result for the electron helicity-flip amplitudes then 
reads 



ImE, 



M 2 
l^ 2 



(30) 



f — 

/ dz (1 — Z) cr 7p (l 



■ 2 ^ > In 



-t 



zt 



2m 



In h - - In ( 1 

m e 2 



(1 — z) 2 m 2 
AE 2 z 



s + M 2 
ImF 5 = — -lmF 4 

S — 



(31) 



Putting these expressions into Eq. I|13|) , we obtain our 
result for the beam normal spin asymmetry, 



B n — 



m e \j —t 



Fx(t) 



An 2 F 2 (t)+rF 2 (t) 



dz (1 — z) <j 1p {w ) { In 



(32) 



-t 



zt 



2m\ In 



A 2 



, 2Ez 1 , , 
In 1 In I 1 

m e 2 



(1 — z) 2 m 2 
4E 2 z 



Before presenting the results of the numerical integra- 
tion, we recall that in the previous work, the cross section 
was assumed to be independent of the photon virtuality 
and roughly constant as function of w 2 , cr 7 . p (?i; 2 , Q 2 ) rj 
<jt = const. In that case, the expression for B n was ob- 
tained [Hfl: 



B„ — 



—tar 



Fi(t) 



4tt 2 F 2 (t) + rF 2 (t) 



In 



(33) 



In Fig. |3 we compare the results for the beam nor- 
mal spin asymmetry of Eq. H32J) with the approximative 
formula of Eq. (|33J) . Comparing the thick and the thin 
lines, we see that the numerical impact of the subleading 
terms in t is very small, and the leading t approximation 
indeed dominates. 



B. Electron helicity-conserving amplitudes 

We next turn to the calculation of the imaginary parts 
of the electron hclicity conserving amplitudes F 2 and F 3 . 
This amounts in calculating the electron phase space in- 
tegral J,. As in the case of J B , we start with the angular 
integration. There are four independent integrals over 
the solid angle in terms of which all other integrals can 
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& 
PQ 




0.2 0.3 0.4 0.5 



Q 2 (GeV 2 ) 



FIG. 2: Beam normal spin asymmetry in the forward regime 
for the lab beam energy Eiab = 3 GeV as function of Q 2 . The 
thick line corresponds to the leading t contribution, while the 
thin line corresponds to the full calculation with subleading 
terms in t. 



be expressed. They are 



h 



dClx 
Q 2 Q 2 



dfli 



i ln(l + Q 2 /m 2 ) 



(34) 



l In (1 + Q 2 /m 2 ) 



w 2 - M 2 + Q 2 
{w 2 — M 2 + Q 2 f 



i In (1 + Q 2 /m 2 ) 

1„ ( A 2 A 



In 

In (1 + Q 2 /m 2 



In 



A 2 



The integrals .Z3 and. I4 do not contain any singularity, 
while one has to be careful with the first two integrals 
in Eq. H34|) . In the limit of small electron mass, one has 
for the first terms in the square brakets (for details, see 

0113) 



h 
h 



2tt 
^tE 2 



In 



E 2 



EEi n m, 



(E- 
AE 2 



Ei) 



(35) 



One can notice that in the limit of zero mass of the 
electron the first two integrals become singular. In the 
case of the beam asymmetry B n , this singularity is can- 
celled by the overall factor m e , thus leading to the domi- 
nant term ~ m e In = \. In the case of helicity conserving 
amplitudes, no such singularity may occur since it does 



not contain the overall factor m e which would cancel the 
chiral singularity. Therefore, the divergent parts of the 
integrals I\ and I2 should cancel exactly in the integral 
I A ■ The integral J A can be written in terms of these four 
integrals as 



{MM™ 2 , 0,0)} 



E 1 dE 1 w 2 — M 2 
I671" 3 27re 2 



<j 1p (w 2 ) 



E 2 + E 2 



v 2 - M 2 2E{E-Ex) 



AE 



E-E x 



h 



1 



E l r 



2K 



2 Ei T 
~E h 



M 2 E 



(36) 



Indeed, in the combination I2 — 2K 2 ^Ii, the depen- 



dence on m P cancels 
h 



2K 2 ^h = 



E ^ EEi "' Q 
Combining all the terms, we obtain: 

{W^^lq 2 )} 



n ln UE-Exf (3?) 



87r 3 e 2 



dz 
1 - z 



a lp {w 2 ) 



IE 

x <j (1 + z 2 ) In — + (1 + z + z 2 ) ln(l - z) 



2E z 



\fs 1 - — ; 



1 



In A 

+ 



M 2 

1 + z 2 



, , M 2 
— z In I 1 z 



2 SP ('^ Z 



1 + z In I i 

2E 



AE 2 



ln- 



1 - z 



V« 1 
-zSp 



M 2 



In 1 



(2E z 
\Ts\~~z 



(38) 



At high energies, the upper limit of the integration 
goes as ^ff- ~ 1 — -j? 2 - — > 1. Consequently, the integrand 
has a rather singular 'behaviour for photon energies justs 
above the pion production threshold. This singularity is 
compensated by the vanishing of the cross section at pion 
production threshold, as it is dictated by unitarity. In the 
high energy approach used here, to ensure the correct 
threshold behaviour we multiply the cross section by a 
function which vanishes smoothly at the threshold and 
goes to 1 at higher energies, such that the high energy 
region is not affected. Here, we will adopt the threshold 
regulator of the following form: 



fthr(z) 




(39) 
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which obeys the properties required above. An introduc- 
tion of such a threshold regulator does, of course, lead 
to a certain model dependence in the low energy part of 
the integral. It is worthwhile to treat low photon energy 
region more accurately, for example, by using some phe- 
nomenological model for the resonant region and match 
it to the high energies contribution, and we leave this 
improvement to an upcoming work. For the moment, 
we concentrate on the high energies-motivated calcula- 
tion. In this framework, the presented approach corre- 
sponds to accounting for the non-resonant background 
which originates from the quark-hadron duality picture. 
Since vanishing of this background function at the pion 
production threshold is required by unitarity, the intro- 
duced threshold function is dictated by physics, though 
not unique. 

Finally, we can write down the analytical part of the 
result for lmF 2 , 



JmFo 



M 2 Q 2 



&ir 2 (s- M 2 ) 



^a 7p (M 2 + (s-M*)(l-z))J^ 



1 - z 



x { (1 + z 2 ) In ^ + (1 + z + z 2 ) ln(l - z) 



2E 



r s i-Ml z 



2 In 1 



M 



r2 



+ - 



In 



A 2 



1 + z 2 



Spl- 



4E 2 



e{k) + e(-fc') -> N(-p) + N(p'). For elastic eN scat- 
tering, the invariant amplitudes possess definite symme- 
try properties between s- and it-channels. These proper- 
ties become more transparent if introducing the explic- 
itly crossing antisymmetric variable v = which is 
positive in the s-channel kinematics and negative in the 
it-channel. At zero momentum transfer this variables re- 
duces to the lab energy of the electron. It can be shown 
that the amplitudes F\^ 2 are odd functions of v, while F 3 
is an even function. As a result, these amplitudes obey 
the dispersion relations of two different forms, 



2v„ 



ReF 1>2 (^) = — V . ,., , 

7T / V U — V 2 



dv' 



Im s F lj2 (v',Q 2 ), 



ReF 3 {v,Q 



■2, l r t i''' 1 "' 



■Im s F 3 (v',Q 2 ), (42) 



where Im s indicates that imaginary parts are calculated 
in the s-channel. We next discuss the convergence of 
these dispersion integrals. As can be easily derived from 
Eq. ©, at very high energies, 



\m. B F 12 {v — > 00, ( 



1 



A'A'^iAAjv 5 



Im s F 3 ^^oo,Q 2 ) -T; 



A'A' ;AAa 



(43) 



where Tyy -x\ N stands for the helicity amplitudes for 
elastic ep-scattering with initial (final) electron (proton) 
helicities indicated. In Regge-regime, the asimptotic en- 
ergy dependence of the latter in the helicity conserving 
channel is given by the pomeron trajectory, 



InxFs 



1 + zln {1 

2E z 
^fsY^MTz 

V S 

(2E z 
-zSp —=- 

M 2 



4E 2 



s-M 2 



ImF? 



In- 



4E 2 
In H ttz 



(40) 
(41) 



T 



A' A^r i A A j\ 



(44) 



with ap ~ 1.08. Inserting these asimptotics in the cor- 
responding dispersion relations for the invariant ampli- 
tudes, we see that the dispersion integrals converge, and 
no subtractions are needed. We therefore can obtain 
clean predictions for the real parts of the invariant ampli- 
tudes from a dispersive calculation which capitalizes on 
the unitarity and analyticity of the invariant amplitudes. 



IV. DISPERSIVE CALCULATION FOR REAL 
PARTS OF F12 3 

In the previous section we obtained imaginary parts 
of the amplitudes Fi, F 2 and F3. Real parts of these 
amplitudes enter the expression of the elastic ep cross 
section. To calculate these, we use dispersion rela- 
tions at fixed t = —Q 2 . Dispersion relations capi- 
talize on the idea of analyticity of the scattering am- 
plitude, which states that the same scattering ampli- 
tude describe the reaction in all three distinct chan- 
nels related by partial CP-transformation, the so-called 
crossing: s-channel e(fc) + N(p) — > e(fc') + A^(p'), u- 
channel e(— k') + N(p) — > e(— k) + N(p'), and f-channel 



V. RESULTS 

We use a numerical routine for the principal value inte- 
grals, and present our results for the real and imaginary 
parts of the elastic ep-scattering amplitudes F 2 and F3 
for 5 different beam energies and as functions of the elas- 
tic momentum transfer Q 2 . To present our results in the 
off-forward kinematics, we need phcnomcnological input 
in order to reproduce the Compton amplitude correctly 
at finite values of t. We use the exponential fit to Comp- 
ton data, 
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FIG. 3: Imaginary (upper panel) and real (lower panel) part FIG. 4: Imaginary (upper panel) and real (lower panel) part 
of the 27-exchange contribution to the amplitude F2. of the 27-exchange contribution to the amplitude F3. 



which provides a good description for B w 7 GeV~2 for 
values of —t up to w 0.8 GeV 2 . Since ^ is related to the 
amplitude squared, while the total cross section to the 
amplitude's imaginary part, we conclude that |l2l | 



a T (t) = a T (t ) x e" 



(46) 



In Figs. 13 and 0] we display results for the amplitudes 
i*2 and F3, respectively, in the forward regime using the 
model described above. The imaginary parts of the two 
amplitudes are related in a simple way, see Ea. (|41|l . and 
therefore are roughly of the same size and of the same 
sign. However, Rei 7 ^ is only about 10% of Rei 7 ^. This 
is due to the different forms of the dispersion relations 
which obey the two amplitudes. In the dispersion inte- 
gral for i<2, V J -jtt-ji. I m -p2 ; the low values of v' ', where 
the denominator is negative, have more impact than in 
that for F3, thus we observe a substantial cancellation in 
this integral. Instead, for F3, the extra power of the en- 
ergy in the numerator, V J -i Imf3 shifts the empha- 
sis on higher values of energy and suppresses low energies, 
therefore the cancellation between the two regions is only 
very moderate. In Fig. we show the relative contri- 
bution of the TPE effects to the form factors weighted 
to the OPE values, 5^lk an d Y 2j = jj^f-- It can 
be seen that the contribution to F2 is fairly small. The 
combination Y 2y is quite small (0.1%) at intermediate e's, 
as well, but grows rapidly with e approaching 1, i.e. at 
higher energies. In this plot, we observe a very steep e 



dependence of the generalized form factors at large e val- 
ues. As one goes to higher momentum transfers, though, 
this behaviour becomes smoother. At very high energies, 
the TPE contributions should decrease, according to the 
high energy asymtotics discussed above. Since we are 
interested in low values of t, such energies would corre- 
spond to the values of e extremely close to 1 and cannot 
be displayed on this plot. 

Finally, we turn to the e + /e _ cross section ratio. In the 
forward kinematics and recalling that in forward regime 
6Gm = 0, we obtain 



R e+/c 



F? + rFi 



Re 



F 1 -^-F 3 -G E SF 2 

1V1 



(47) 



The numerical results are shown in Fig. |SJ In the for- 
ward regime, the positron cross section is expected to 
exceed that for electrons by a few procents: from 1.4% 
for Ei ab = 3 GeV to 2.8% for E lab = 45 GeV. If going 
to even higher beam energies, the ratio drops due to the 
high energy asymptotics of the invariant amplitudes. At 
the lab energies as high as several tens GeV, the disper- 
sive contributions presented here are expected to be the 
dominant effect, while the elastic contribution (nucleon 
intermediate state in the blob in Fig. ^) should be rela- 
tively small. So, the above results should give a correct 
estimate of the TPE contribution to the elastic cross sec- 
tion in forward regiime. 
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We presented a dispersive calculation for the real parts 
of the electron helicity-conserving amplitudes -Fi,2,3 for 
elastic e ± p-scattering in forward regime. Their imagi- 
nary parts in this regime can be related to the DIS struc- 
ture functions. We used the phenomenological w 2 and 
Q 2 dependence of the structure functions. The imaginary 
parts of the ep scattering amplitudes are IR finite, and we 
demonstrated that their high energy asymptotics ensure 
the convergence of the unsubtracted dispersive integrals 
for the corresponding imaginary parts. Due to differrent 
forms of the dispersion relations for the amplitudes Fi 
and i*3, the resulting real parts have very different val- 
ues: the dispersive effects contribute only about 0.05% of 
the low Q 2 -value of F2, while they account for a half to 
few per cents in the case of Y27 = jj ^Q Fa ■ We predict the 
ratio of the e + p to e~p cross sections to be larger than 1, 
and the predicted values for that ratio range from 1.4% 
at Ei ab = 3 GeV to 2.8% at E lab = 45 GeV. 
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FIG. 5: Real part of 6F2 normalized to the OPE value of F2 
at three different values of Q 2 and as function of e (upper 
panel). The same for the 27-amplitude Y27 (see text for the 
definition) is displayed in the lower panel. 
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FIG. 6: e + /e asymmetry in the forward regime for five dif- 
ferent values of the lab beam energy as function of Q 2 . 



